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Abstract. ESO 338-IG04, also known as Tololo 1924-416, is a well known luminous blue compact galaxy in the 
local universe. Images obtained with the Hubble Space Telescope (HST) have shown that the central starburst 
region is composed of numerous bright point sources - young star clusters, surrounded by a population of old and 
intermediate age globular clusters. In this paper we use ultra-violet (UV) and optical HST photometry in five 
bands, and an extensive set of spectral evolutionary synthesis scenarios to investigate the age and masses of 124 
star clusters. The very small reddening makes ESO 338-IG04 an excellent laboratory for studying the formation 
of such objects. We find that a careful treatment of the nebular emission component is crucial when modelling the 
broad-band colours of young starburst regions. We have used the star clusters to trace the temporal and spatial 
evolution of the starburst, and to put constraints on the star formation activity over a cosmological time-scale. 
The present starburst has been active for about 40 Myr and shows evidence for propagating star formation and 
structures triggered by galactic winds. A standard Salpeter initial mass function (IMF) extending up to 120 Mq 
provides the best fit to the data, although a flatter IMF cannot be excluded. The compact star clusters provide 
30-40% of the UV luminosity and star formation activity. We find no evidence for dust obscuration even among 
the youngest (< 1 Myr) clusters, and we propose that this may be related to a short time-scale for destruction 
of dusty molecular clouds. Over a longer time-scale, we find evidence for previous cluster formation epochs - 
notably one a couple of Gyr ago. The fraction of the galaxy's stellar mass contained in compact star clusters is 
found to be several percent, which is an unusually high value. The intermediate age clusters show a flattened 
space distribution which agrees with the isophotal shape of the galaxy, whereas the oldest clusters seem to have 
a spherical distribution indicating that they formed prior to the rest of the galaxy. 

Key words. Galaxies: starburst - galaxies: star clusters - galaxies: compact - galaxies: individual: ESO 338-IG04 
(= Tol 1924-416) 



1. Introduction 

According to the hierarchical paradigm for structure for- 
mation, small structures form first and massive galaxies 
are successively built up by mergers of smaller units. Blue 
Compact Galaxies (BCGs) are galaxies with lower mass 
and metallicity than typical L* galaxies. Their spectra 
often show a blue continuum and strong emission lines, 
indicative of low chemical abundances and active star for- 
mation. These properties once led to the idea that BCGs 
were young systems presently forming their first genera- 
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tion of stars (Searle & Sargent Il972[l . Today it is evident 
that the majority of BCGs are not young (see Kunth & 
Ostlin 2000 for a review) , but rather experiencing a phase 
of increased star formation activity, ranging from moder- 
ate enhancements (Schulte-Ladbeck et al. 2001, Hopkins 
et al. l2002|l to true starbursts, with short gas consumption 
time-scales (Ostlin et al. 120011) . Understanding the reasons 
for this lies at the heart of the quest to understand the role 
of BCGs in the cosmic galaxy evolution framework. 

Observations of e.g. the Hubble Deep Field (Williams 
et al. 11996(1 . show that luminous compact blue galaxies 
are common at higher redshifts (Guzman et al. ^997). 
Moreover, Lya emitters at redshift z — 2.4 have lower 
continuum luminosity and smaller sizes than expected for 
young L* galaxies (Pascarelle et al. :1998). These may 
be dwarf galaxies, about to build up larger structures 
by mergers. Several investigations (e.g. Lilly et al. [1996 
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Madau et al. fTMrJj) show that the cosmic star formation 
rate (SFR) was much higher at z > 1, and the same ap- 
plies to the merger frequency although mergers appear to 
have played an important role also at z < 1 (e.g. Le Fevre 
et al. I2000[l . Based on kinematics, morphology and pho- 
tometry, Ostlin et al. (200! ) suggest that mergers may be 
the main mechanism responsible for starbursts in lumi- 
nous BCGs (LBCGs) in the local universe. This is also 
supported by the enigmatic under-abundance of the ISM 
metallicity compared to the stellar population found by 
Bcrgvall & Ostlin 1|2002[1 . Hence, LBCGs offer an oppor- 
tunity to study hierarchical buildup of galaxies in detail. 

Globular clusters (GCs) include the oldest known stel- 
lar systems in the local universe. Old Galactic halo GCs 
are around 15 Gyr old and typically have metallicities 
[Fe/H]= -1.5 (Ashman & Zepf [TU^jl . High resolution 
imaging, notably with the Hubble Space Telescope (HST) 
of merging galaxies in the local universe have shown that 
luminous compact star clusters (often referred to as su- 
per star clusters) with properties consistent with young 
GCs are formed in such events (see Whitmore 120021 for 
a recent review). This can explain the existence of in- 
termediate age (a few Gyr) metal-rich GC populations 
found in many local elliptical galaxies (see Ashman & 
Zepf 1 199*51 for a review). Hence, GC formation was not 
confined to the high z universe, but is still going on but 
at a slower pace. Also LBCGs are sites of young compact 
star clusters (Ostlin 2000). However, the ISM abundances 
([0/H]« —1) are lower than those in giant mergers, and 
similar (if [Fc/0]~ —0.5, as is common in metal-poor envi- 
ronments, see e.g. Kunth and Ostlin 2000J) to the metallic- 
ity of old GCs. The first generation of GCs in the universe 
may have formed in the processes of merging sub-galactic 
clumps, similar to what is going on in some local LBCGs. 

ESO 338-IG04 is a well known luminous (My = -19.3) 
BCG, first studied in detail by Bergvall (1985 ). Images ob- 
tained with the Wide Field Planetary Camera 2 (WFPC2) 
on board the HST found numerous compact star clusters 
associated with ESO 338-IG04 (Ostlin et al. ITM81 here- 
after OBR98). In addition UV observations with the Faint 
Object Camera (FOC) on board HST, had already shown 
that the centre of ESO 338-IG04 hosted a population of 
UV-bright star clusters, just like many other starbursts 
(Meurer et al. H995|l . Further analysis showed that the star 
clusters in ESO 338-IG04 spanned ages from a few Myr to 
more than 10 Gyr and that they had just the properties 
expected for young and old GCs (OBR98). The galaxy is a 
true starburst in the sense that the time-scales for gas con- 
sumption and build-up of the observed stellar mass with 
the current SFR both are ~ 1 Gyr, i.e. much shorter than 
the age of the universe (Ostlin et al. 2001). 

In the present paper, we will make a refined analysis 
of 124 star clusters in ESO338-IG04 using additional UV 
data and a new set of models. We use the star clusters and 
their age- and space-distributions to track the evolution of 
the star formation in the galaxy over a cosmological time 
scale. 



2. Observational data 

ESO 338-IG04 was observed with HST/WFPC2 with the 
galaxy centred on the Planetary Camera (PC) aperture 
(OBR98). The photometry in the F336W, F439W, F555W 
and F814W bands (hereafter referred to as: u, b, v, i) is de- 
scribed in OBR98. In addition we use photometry in the 
F218W filter which is very valuable for constraining the 
physical parameters and the extinction of young objects. 
The F218W photometry (hereafter x) has been performed 
in the same way as for the other passbands, and has 
been corrected for the same amount of Galactic extinction: 
E(B - V) G = 0.09 (Burstcin and Heiles 1982). We have 
also corrected the F218W data for a "UV-contamination" 
of -0.1276 magnitudes (Biretta et al. E002jl . All photom- 
etry in this paper is in the VEGAMAG system (Biretta 
et al. l2TjU5|l 

Figured] shows an RGB true colour composite of the i, 
v and b images. The v filter (F555W) has significant trans- 
mission at H/3, [O in] A4959, 5007 and even Ha. Consequently, 
several nebulous emission-line regions are seen in green 
in Fig. n It is obvious that this galaxy presents a large 
number of compact star cluster candidates. In OBR98 we 
showed that the vast majority of these are indeed star 
clusters, physically associated with ESO338-IG04 (i.e. not 
foreground stars, background galaxies or individual super- 
giants in ESO338-IG04). Their sizes are compatible with 
globular clusters, although a few of the younger sources are 
slightly more extended and may represent unbound sys- 
tems. In fig. [2] we show an RGB composite of the central 
region using the i, v and x filters. The variety of colours ev- 
ident from Figs.n an d0 a central population with bright 
blue clusters surrounded by a population of fainter red 
ones, directly suggests that the clusters have different age 
or reddening. 

Figure|3]shows the F218W magnitude (m x ) vs the pho- 
tometric uncertainties (cr x ), and vs M v , the absolute v 
magnitudes for the same distance as used in OBR98 (37.5 
Mpc, or m-M = 32.87) 1 . From the catalogue in OBR98, a 
total of 78 objects with m x < 22 were detected. However, 
cosmic ray contamination is worse (due to longer expo- 
sure time, 1800s split into two sub-exposures) than in the 
other filters, and some spurious detections had to be re- 
moved by hand. In total, 49 objects have a x < 0.5 and 
21 objects have a x < 0.2. This is comparable to Meurer 
et al. JT995) who detected 14 compact sources with the 
pre-refurbishment HST/FOC and the F220W filter hav- 
ing uncertainties <tf220W < 0.17. The photometric agree- 
ment is fair when comparing sources for which Meurer et 
al. I|1995|l determined the flux from radial profiles. From 
Fig. |3 it is evident that the sources detected in F218W 
have very blue x — v colours, with a median around —2, 
a combined effect of low age, metallicity and extinction of 
the UV detected clusters. 

Figurc^Jshows the x — uvsb — v two-colour diagram for 
objects with accurate photometry, together with two typ- 

1 The distance is based on a Hubble parameter of Ho = 75 
km/s/Mpc, which is assumed throughout this paper. 
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Fig. 1. See separate jpg-file: 3277.fl.jpg True colour image composite of the central 30" x 30" (corresponding 
to 5.45 x 5.45 kpc at the adopted distance 37.5 Mpc) region of the Planetary Camera. F814W (i) is coded in red, 
F555W (v) in green and F439W(&) in blue. North is up-left (65 degrees measured from vertical towards left) and east 
is down- left. The bright source with diffraction spikes ~ 5" from the centre towards the upper right is a Galactic 
foreground star. 



Fig. 2. See separate jpg-file: 3277. f2.jpg True colour 
RGB composite of the central 11.6" x 11.6" (2x2 kpc) re- 
gion with cuts chosen to show details of the central bright 
regions: F814W (i, red), F555W (v, green) and F218W(x, 
blue). The greyscale inset shows the F218W image in of 
the same field. 
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Fig. 3. Photometry for clusters detected in the F218W 
(x) band. The bottom panel shows the magnitude- 
uncertainty diagram. Top panel shows the apparent x 
magnitude vs the absolute v band magnitude. Lines of 
constant colour are shown as well: x — v = (dotted), —2 
(dashed), and —4 (dash-dot). 

ical models used (see below for details on these). During 
the first 10 Myr the x — u colour is nearly constant making 
it an efficient probe of the reddening for young objects. On 
the other hand, b — v evolves strongly, mainly due to the 
rapid decrease in the contribution from emission lines to 
the v (F555W) filter. That the objects scatter close to the 
tracks demonstrate that the extinction is very small. It is 
also clear that it is only for objects younger than a few 
times 10 Myr that we have reliable F218W data. 

In addition to the clusters in Fig. 2] which are found 
in the unvignetted field of view (FOV) of the PC, a blue 
tail of the galaxy (see Bergvall & Ostlin I2002H extends 
into the WF4 aperture (downwards in the orientation of 
Fig.Q, where another dozen compact objects are found. 



This is near the border between the chips, and we do not 
regard the photometry of these objects reliable. Some ob- 
jects appear extended, others resemble compact star clus- 
ters, which is confirmed by our new deep STIS data (Ostlin 
et al. in prep.). The STIS images (FOV 9.3 x 9.3 kpc) also 
show a small number of bright clusters lying outside the 
upper and right boundaries of the PC FOV. Bergvall and 
Ostlin trace the luminosity profile of ESO 338-04 to radii 
greater than 10 kpc, and in OBR98 it was shown that the 
radial cluster distribution follows that of the visual light. 
From this we estimate that the limited area of the PC re- 
sults in a cluster incompleteness, at the ~ 11% level, in 
addition to the photometric incompleteness. 

In Table |31 we present the integrated starburst lumi- 
nosity in different passbands and the fraction of light at- 
tributable to compact star clusters. 

3. Spectral evolutionary synthesis models 

Due to the uncertainties still prevailing the field of spectral 
evolutionary synthesis, all results should ideally be cross- 
checked with independent models. To ensure that the age 
and mass distributions derived for our star cluster sample 
do not critically depend on the predictions from any single 
code, the analysis has been carried out with two models 
in parallel: PEGASE.2 (Fioc & Rocca- Volmerange HSGfy 
and the model by Zackrisson et al. 1)20011 hereafter Z2001). 
These particular models were chosen since they contain a 
nebular component and because of the ease with which 
non standard filters may be included. 

PEGASE.2 is based on the method of isochrone syn- 
thesis, employs stellar atmospheres by Lejeune (|1998|) and 
Clegg & Middlemass (|1987JI . together with stellar evolu- 
tionary tracks mainly from the Padova group. Using a 
simple prescription for the nebular component, it predicts 
both continuum and line emission of the ionised gas. 

Z2001 is based on the method of isomass synthesis, em- 
ploys the same set of stellar atmospheres as PEGASE.2 
but stellar evolutionary tracks mainly from the Geneva 
group. Pre-main sequence evolution and a stochastic treat- 
ment of horizontal branch morphologies at low metallici- 
ties are also included. Gas continuum and emission lines 
are predicted using the photoionisation code Cloudy ver- 
sion 90.05 (Ferland 1996|. For each time step, the spec- 
tral energy distribution of the stellar population is fed into 
Cloudy, making the computation time consuming, but also 
giving a more realistic nebular emission component. 

The use of pre-main sequence evolution and and a 
more sophisticated nebular component in Z2001 makes 
this model more suited for the current analysis, where 
many clusters are believed to be young. As will be demon- 
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Fig. 4. Colour-colour diagram including the F218W data 
(denoted 'x'). The median colour uncertainty for objects 
included in the plot is shown as a cross in the lower right, 
whereas the largest uncertainty associated with any ob- 
ject is indicated by the surrounding rectangle. The ar- 
row indicates how much and in what direction the ob- 
jects would be repositioned after correction for a hypo- 
thetical extinction of E(B — V) = 0.25. The predictions 
for two of the used spectral evolutionary synthesis mod- 
els (Salpeter IMF, M = 0.08 to 120M Q , instantaneous 
burst, Z gas = 0.002, Z sta rs = 0.001) are shown as well: 
Zackrisson et al (2001, solid line), and PEGASE.2 (dashed 
line). Markers indicate ages of 0.5 Myr (circles), 10 Myr 
(triangles), 1 Gyr (squares) to 15 Gyr (pentagrams). 



strated (see Sect. 4.1, 5.1) the predictions by Z2001 are 
also in better agreement with the observations. A large 
grid of evolutionary sequences was therefore generated us- 
ing Z2001, and a smaller set, intended to use as a con- 
sistency check, with PEGASE.2. The details of these se- 
quences are given in tables 1 and 2. In total, 133 param- 
eter configurations were processed with Z2001, 10 with 
PEGASE.2. The two models used represent an improve- 
ment over the one used in OBR98 in terms of the nebular 
component, time resolution at young ages, updated stel- 
lar evolutionary tracks and stellar atmospheres at a wider 
range of metallicities. In addition, the Z2001 grid repre- 
sents an increase in the size of the parameter space by 
more than a factor of 30 compared to OBR98. 

Magnitudes in the HST Vegamag system were gen- 
erated from model spectra using the HST/WFPC2 fil- 
ter+CCD response curves, and a theoretical Vega model. 
Since the full transmission profiles were used, red leaks 
are included in the calculation and do not result in any 
degenercies. The stellar initial mass function (IMF) was 
parametrised as a power-law: dN/dM oc M~ a , with spec- 
tral index a. For a Salpeter <fT§5l)|l IMF, a = 2.35. The 
Scalo (1998, hereafter Scalo98) IMF used in some of the 
PEGASE.2 evolutionary sequences was parametrised by a 
piecewise continuous 3 segment power-law. 



Table 1. The grid of Z2001 (Zackrisson et al. I^OlJTjl 

evolutionary sequences. The grid consists of all possi- 
ble combinations of the parameter values listed below. 
To furthermore ensure the robustness of the final re- 
sults, a smaller number of evolutionary sequences was 
also calculated with variations like: constant star forma- 
tion during 5 Myr and 100 Myr, M up = 40, 60, 8OM , 
10 6 A/ Q available for star formation, n(H) = 100 cm -3 , 
filling factor 1.0, Z stars = 0.004, 0.008, 0.020, 0.040 and 
Z gas = 0.001, 0.004, 0.008, 0.020, 0.040. Neither the vari- 
ations in hydrogen density, filling factor nor stellar mass 
were seen to have any serious impact on the final ages and 
masses derived. 



IMF, a 

Mlow 

M U p 
SFH 

■^stars 
■^gas 

ars 

n(H) 

Filling factor 
Covering factor 
Nebular emission 



1.35, 2.35, 2.85 

0.08, 2 M Q 

20, 120 M Q 

cle7, e3e6, Inst 

0.001 

0.002 

10 5 M Q 

10 cm~ 3 

0.1 

0.5, 1.0 
Yes, No 



IMF: dN/dM oc M~ a 

SFH=Star formation history: cle7 = Constant SFR during 
10 Myr, e3e6 = Exponentially declining SFR with e-folding 
decay rate of 3 Myr, Inst = Instantaneous burst. 



Table 2. List of used PEGASE.2 evolutionary sequences, 
all assuming ^stars — ^gas — 0.002. All evolutinary se- 
quences were computed both with and without inclusion 
of nebular emission. 



IMF 


M low (M ) 


M up (M Q ) 


SFH 


Salpeter 


0.1 


120 


Inst 


Salpeter 


0.1 


20 


Inst 


Salpeter 


0.1 


120 


e3e6 


Salpeter 


1 


120 


Inst 


Scalo 1998 


0.1 


120 


Inst 



4. Modelling the ages of star clusters 

The object magnitudes (corrected for aperture effects, 
Galactic extinction, and UV-contamination) were con- 
verted into fluxes and assigned weights inversely propor- 
tional to the square of the uncertainties. Since the un- 
certainties are based on the photon statistics, and there 
are other error sources present as well, e.g. zero-point and 
flat field errors, a maximum weight corresponding to an 
uncertainty of 0.03 magnitudes was used. The contribu- 
tion from emission lines is significantly larger in the v fil- 
ter (F555W). Since the nebular component is sensitive to 
the assumed ISM parameters, and because of the diffi- 
culty in isolating all the diffuse nebular emission associ- 
ated with a particular star cluster, we increased the formal 
uncertainties with a factor of 2 in this filter. The predic- 
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Table 3. Integrated starburst luminosities. All values 
have been corrected for foreground reddening by our 
Galaxy. The first column gives the name of the filter, 
and the second column the effective wavelength of the 
filter+CCD combination (Biretta et al. 2002). The third 
column gives the integrated starburst magnitude (in the 
VEGAMAG system) within a radius of 5.9"(1.075 kpc) 
centred on the brightest cluster. The fourth column gives 
the same quantity in f v units. The fifth column gives the 
magnitude of the starburst region after all detected com- 
pact clusters have been subtracted. The last column gives 
fractional contribution from compact star clusters to the 
integrated starburst luminosity. 



Filter 




Starburst 




Diffuse 


Lci/Lsb 






mag 


u 




mag 






A 


vega 


W/m 2 /Hz 


vega 




x (F218W)" 


2210 


12.61 


6.6 10" 


29 


12.97 


0.29 


u (F336W) a 


3348 


13.35 


5.7 10" 


-29 


13.66 


0.25 


6 (F439W) 


4316 


14.58 


6.3 10" 


-29 


14.83 


0.21 


v (F555W) 


5465 


14.08 


8.7 10" 


-29 


14.31 


0.19 


i (F814W) 


8040 


14.07 


5.9 1(T 


-29 


14.37 


0.15 



Note: a - corrected for UV-contamination 0.12 mag in F228W 
and 0.01 mag in F336W 



1.6 
1.4 
1.2 

1 

0.8 



> 0.6 



0.4 
0.2 



-0.2 
-0.4 




-0.5 




v-i 



0.5 



Fig. 5. Colour-colour diagram (v — i vs b — v), see caption 
of Fig. 0] for more details. The full drawn line show the 
prediction form the Z2001 model with Salpeter IMF, M = 
0.08 to 12OM0, instantaneous burst, covering factor = 1.0, 
Z gas = 0.002, Z s tars = 0.001. The dashed line shows the 
prediction from the publicaly available PEGASE.2 model 
with similar model parameters (instantaneous burst with 



0.08-120M Q Salpeter IMF, Z gas = 0.002, Z s , 



0.002). 



tions from the spectral synthesis models, having 51 time 
steps for the Z2001 model, is interpolated linearly onto a 
finer grid with an age resolution of 1 Myr up to 3 Gyr 
and 100 Myr from 3 to 15 Gyr. Finally, we find the best 



fitting age for each cluster by minimising the weighted 
RMS deviation. When allowing E(B ~ V) to vary, we used 
E(B - V) = 0.0,0.05,0.10,0.15,0.20 and 0.25 and deter- 
mined the best fitting value for each cluster. 

Masses were determined by scaling the extinction cor- 
rected v flux to the M/L v predicted by the model. These 
take into account the total mass converted into stars 
and thus accounts also for stellar remnants and gas re- 
turned in the process of stellar evolution. The "true" stel- 
lar mass (excluding returned gas but including remnants) 
will therefore be overestimated at high ages, but since we 
are interested in the star formation history, we prefer us- 
ing the total mass initially present in each cluster. At 15 
Gyr the overestimate in M/L v will amount to a factor of 
-~ 1.5 for 0.08 - 120M Q Salpeter IMFs (for any Z stars ), 
~ 3 for a = 1.35, and ~ 1.1 for a = 2.85. 

In Fig. [S] we show the v — i vs b — v diagram for clus- 
ters with errors smaller than 0.15 mag in both colours to- 
gether with a predictions from the Z2001 and PEGASE.2 
codes for a standard set of parameters. Only Z2001 is suc- 
cessful in reproducing the youngest clusters, emphasizing 
the need for a careful treatment of the nebular contribu- 
tion. For ages larger than a few times 10 Myr, both mod- 
els approximately perform equally well. However, around 
a few 10 Myrs the PEGASE.2 model presents loops in 
the colour-colour diagram (Fig. 03) which the Z2001 model 
does not. This different behaviour is due to the use of dif- 
ferent evolutionary tracks (see Chariot et al. 119961 for a 
discussion). 

4.1. Internal reddening 

The effects of age and reddening are in general difficult 
to disentangle, although the reddening and aging vectors 
are not absolutely parallel. In particular, the reddening 
vector for young objects is close to orthogonal to the age- 
metallicity sequence (Fig. EJ) and, as has also been shown 
by other investigators (e.g. Mas-Hesse and Kunth, 1999), 
the addition of 2000A data helps in deriving the extinc- 
tion for very young objects, since the UV-spectral slope is 
nearly constant at low ages (see Fig. 0J. 

In ESO338-IG04 we have several reasons to suspect 
that the internal reddening is very small. Groundbased 
studies of the Balmer emission line decrement (Bergvall 
1986, lye et al. 1987, Calzetti et al. I1HH4| Raimann et al. 
2000), as well as studies of the UV continuum (Calzetti 
et al. IT5M1 Meurer et al. Buat et al. IMHSI of the 

central star forming region have given values of E(B — V) 
in the range from to 0.05. Further support from this 
conclusion comes from a comparison of the far IR to UV 
luminosity (Meurer et al. 119991 Buat et al. E002) and the 
colours of the bluest star clusters themselves (this paper) . 
Hence, we know that the average extinction is very small, 
making ESO338-IG04 an ideal target to study the forma- 
tion and evolution of star clusters. However, star forma- 
tion is normally associated with dusty molecular clouds, 
and it is likely that the reddening varies from cluster to 
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cluster. Especially we may suspect the existence of embed- 
ded very young sources, that, however, will have quite low 
probability of being picked out from this optically selected 
sample. Groundbased long-slit spectra, obtained with the 
ESO NTT shows a Ha/H/3 value that varies along the 
slit, corresponding to a varying E(B — V) in the range 
0.0 to 0.20, although the luminosity weighted average is 
close to zero. To adress the internal extinction we have 
adoped two approaches: i) Allowing E(B — V) to be a 
free parameter in the fit, but constraining it to the range 
E(B — V) £ [0,0.25]. ii) Adopting a single extinction for 
all clusters in the range E(B - V) G [0, 0.25]. 

In general we use a Galactic extinction law (Seaton 
I1979f) . which has an intermediate behaviour in the UV 
compared to the LMC and Calzetti laws commonly used 
for starbursts (Calzetti 1994). In the optical, the difference 
is small. Tests perfomed with the Calzetti law give close 
to identical results regarding age and best fitting E(B — 
V). Hence, the results presented are not sensitive to the 
choice of extinction law within the range of low extinctions 
considered. However, for modelling the UV emission of 
clusters in more reddened galaxies, the extinction law will 
be an additional source of uncertainty. 

4.2. Comments on the parameter choices 

In most model runs, we have kept the parameter choice 
fixed, but in some cases, we have allowed some parameters 
to be free, as e.g. the covering factor or the IMF slope. 
There is accumulating evidence for an IMF slope very 
close to the Salpeter value for stars more massive than 
1M Q (Massey & Hunter [TUM Kroupa I2UU21 Kennicutt 
1998, Leitherer 1998). At lower masses, the mass function 
seems to turn-over with the result that the IMF can be 
approximated by a multi-segment power-law (e.g. Miller- 
ScaloEHS Scalo 1998). Since the main sequence life time 
of a IMq star is on the order of 10 Gyr, the shape of the 
lower mass IMF will have very little influence on colours 
for any age considered in this paper. The main effect lies 
in a different mass normalisation leading to lower effective 
M/L values for the Scalo98 fScalo. HMs) IMF for ages 
greater than ~ 300 Myr, compared to a Salpeter IMF. 

The choice of IMF slope (a), upper mass limit and 
covering factor has important effects on the colours at low 
ages. At high ages, the most important effect is on the 
mass to light (M/L) values, introducing uncertainties in 
the masses estimated from photometry. 

4.2.1. Covering Factor 

In our Z2001 model runs we have used two different gas 
covering factors in the calculation of the gas spectrum: 1.0 
and 0.5, describing what fraction of ionising photons along 
the line of sight that will be absorbed. This influences the 
strength of emission lines relative to the continuum, and 
hence the predicted flux, notably in the v filter. A covering 
factor of 0, corresponding to no nebular emission contribu- 



tion at all, was also considered, but could not fit the data 
for young/blue objects. The choise of 0.5 is somewhat ar- 
bitrary, but was included to have an intermediate value to 
get an estimate of the importance of this parameter. 

Figure shows that a fraction of the ionised gas emis- 
sion is found outside the bright clusters, but this is ex- 
pected since even in the UV, young clusters don't pro- 
duce more than 30-40% of the starburst luminosity. The 
ionising continuum from a young cluster will produce sig- 
nificant nebular emission, that may have a different spatial 
distribution compared to the clusters. Narrow-band [O in] 
images obtained with HST/STIS (Ostlin et al. in prep.) 
show that the peaks in the [Oiii] surface brightness in- 
deed coincide with blue UV-luminous sources. However, 
the [O in] is more diffuse than the stellar continuum, and 
our apertures could miss a fraction of the gas emission. 
Due to the ionisation structure, a larger fraction of the 
[O in] emitting zone will be included as compared to the 
[Oil] and Hll zones. Since our apparently youngest clus- 
ters still have as blue (v — i) and red (b — v) colours (where 
the v band excess mainly comes from [O in] line emission) 
as the earliest model steps, the net contribution from gas 
outside the aperture must be small. 

Maiz-Appellaniz ( 2001 ) studied the spatial structure of 
gas emission in a sample of young star clusters, and found 
that ~2 Myr old clusters have small (~10 pc) shells around 
them, whereas ~4 Myr clusters have larger shells, and for 
older clusters only diffuse Ha emission is detectable. This 
indicates that with photometric aperture diameters of ~50 
pc, we would not miss any emission during the first few 
Myr, and the covering factor should be close to unity. At 
later times (5-10 Myr), we may miss a significant fraction 
of the emission, which can be approximated by a smaller 
covering factor. At still larger ages the clusters will be too 
cool to ionise [O in] and the ionisation structure will be of 
little concern. Moreover, the relative gas contribution will 
be almost negligible in all passbands used here. 

4.2.2. Metallicity 

The nebular oxygen abundance in ESO338-IG04, as deter- 
mined using the temperature sensitive [Oiii]a4363 line, has 
been shown to be on the order of 10-14% solar (Bergvall 
1984, and Masegosa et al. 1994). Hence, in general, we 
adopt Z gas = 0.002, (Z Q = 0.0188). Other values (see 
Table P) are also included to test how the assumed nebu- 
lar abundances affect the results. It is comforting that the 
best fit is provided by Z gas = 0.002. 

For the stars, we have very limited information on the 
metallicity. In general Fe, which is the dominant source 
of stellar opacity, is less abundant (compared to solar) 
than O in metal-poor galaxies (Kunth and Ostlin 2000). 
Also, in a closed box model, the stars should on aver- 
age have lower metallicity than the gas. Hence, the choice 
^stars = 0.001 seems reasonable. However, there are in- 
dications from near-IR surface photometry that the un- 
derlying old host galaxy may be more metal rich than 
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the nebular gas (Bergvall and Ostlin 2002). This intrigu- 
ing result can be interpreted as metal-poor gas being ac- 
creted onto a more metal-rich host galaxy in a merger. 
Whereas young objects should have a metallicity close to 
that of the gas, it cannot be excluded that older clus- 
ters are significantly more metal-rich. However, metal- 
rich elliptical galaxies have in general metal-poor (typi- 
cally [Fe/H] < -1, Brodie and Huchra lTWT|) old globu- 
lar cluster systems. The commonly found metal-rich sub- 
populations are in general younger and still less metal-rich 
than the host galaxy (Ashman & Zepf 1998). 

4.2.3. Star formation time-scale 

The assumed time-scale for star formation tsf, i.e. how 
long star formation is active in a particular event, impacts 
the model predictions for ages comparable to tsf- As long 
as star formation is active, young stars will make colours 
blue and modelled ages tend to increase with the assumed 
tsf- This is a major source of uncertainty when modelling 
the integrated light from galaxies. In these compact (r ~ 
a few pc) star clusters we know that tsf must be short, 
which is a main reason why ages of such are less degener- 
ate than those for unresolved populations. Of course, the 
common assumption of a truly instantaneous formation 
of a massive star cluster is unrealistic. Another question 
is whether low- and high-mass stars form simultaneously, 
or with some time lag inbetween. When massive stars ap- 
pear, they ionise their surroundings and inhibit further 
star formation. Hence, low mass stars must form prior to 
or simultaneously with the massive stars. McKee and Tan 
(2002) argues that the formation time scale for individual 
stars is short, of the order of 10 5 years, permitting the 
coeval formation of low and high mass stars in a cluster. 
Massey and Hunter (1998) found that the formation of 
R136, the central cluster in 30 Doradus, might have been 
extended over ~4 Myrs, with low mass stars forming prior 
to massive ones. Moreover the IMF slope is a — 2.3 — 2.4 
over the studied range 2.8-120 M@, i.e. identical to the 
Salpeter (1955) value. Similarly, the central star cluster 
in the Galactic giant Hn region NGC3603 shows an age 
spread of a few Myr (Pandey et al. 2000). 

Our three main SFHs considered are instantaneous 
burst, a constant SFR for 10 Myr after which the SFR 
drops to zero, and an exponentially declining SFR with 
an e-folding time of 3 Myr. A smaller number of other pa- 
rameter options were tested as well (Table A potential 
drawback with the exponential model is that the SFR only 
goes asymptotically to zero, and hence remains significant 
for quite a long time. Although little of what is derived 
in this paper will be substantially changed whether one 
adopts either of these three timescales, we will show be- 
low that comparison with the global Ha output tend to 
favour the shorter one. An instantaneous burst also gives 
on average lower RMS residuals. 

In all our evolutionary sequences, low- and high- mass 
stars form simultaneously. If in reality there is a time lag 
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Fig. 6. Some results from the Monte-Carlo simulations. 
Top row: Results from a typical evolutionary sequence 
(Salpeter IMF, M = 0.08 to 12OM , 10 Myr burst 
with constant SFR, covering factor 1.0, Z gas — 0.002, 
^stars = 0.001). On the left we show the age distribution, 
and on the right the E(B — V) distribution. This age and 
E(B — V) distribution was used to create a fake catalogue 
where the objects match the evolutionary sequence per- 
fectly in all filters. To this fake catalogue, random errors 
were added according to the observed uncertainties, and 
then used as input to the fitting program. Second row: 
Sample output from this experiment when E(B — V) was 
a free parameter. In the E(B — V) histogram, the "in- 
put" E(B — V) distribution of the top panel is shown as a 
thick unshaded histogram. Third row: Output E(B — V) 
distributions from sample Monte-Carlo simulations where 
the input fake catalogue was created with fixed values of 
E(B — V). The input E(B — V) value is given as a num- 
ber in each sub-panel, and the E(B — V) distribution of 
the top panel is shown as a thick unshaded histogram for 
reference. 



of a few Myr in the formation of the most massive stars, 
this will influence our results for ages on the order the time 
lag. For a time lag of 3 Myr our model would behave as 
follows: Looking at Fig.[S]for reference, the time-lag tracks 
would start near v — i = — 0.1, b — v = and climb along 
the simultaneous tracks up to v — i — —1.2, b — v = 1.2, a 
point it would reach when the most massive stars appear 
after 3 Myr. Then the spectrum would essentially follow 
the same track as for the simultaneous case model. Hence 
the only effect would be that we would confuse ages during 
the first 3 Myr. We will discuss the possible implications 
of time-lag on our results when relevant. 
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4.3. Monte-Carlo simulations 

In order to assess the internal accuracy of the fitted ages 
and extinctions, we constructed catalogues of artificial 
clusters with the same age and luminosity distribution as 
indicated by the model fits and with colours that per- 
fectly fit the models. To these ideal clusters we added ob- 
servational errors in accordance with the observed error 
distribution. Then this fake catalogue was fitted to the 
models. The experiment was repeated 100 times to beat 
down the statistical noise and from the dispersion of the 
fitted ages we could for each object assess the uncertainty 
in the derived age, for the assumed parameter configura- 
tion in question. The top two panels in Fig. [fj] show the 
input (top) and output (second row) age and E(B — V) 
distributions for a sample Monte-Carlo experiment. The 
tendency to slightly overestimate the output extinction is 
statistically significant, but is only a few hundredths of a 
magnitude on average. 

We also created fake catalogues with fixed input E(B — 
V) values. Specifically, we wanted to test whether the 
E(B — V) range found from the free fit could be an arti- 
fact produced from an (in reality) single reddening value 
smeared out by observational errors. The bottom panel 
in Fig. compares the reddening found from the free 
fit method to simulations with fixed input E(B — V). 
Obviously, no single reddening value can be found which 
can reproduce the distribution found when letting E(B — 
V) be a free parameter, although E[B — V) = comes 
quite close. Hence, we conclude that there is a real spread 
in E(B — V), but letting it be a free parameter on average 
may lead to a slight overestimate of the true E(B — V). 

Since the colour and luminosity of a cluster changes 
with time, the uncertainty in the derived ages is not time- 
independent. We have used the Monte-Carlo simulations 
to estimate the internal uncertainties on the derived ages 
resulting from observational errors. The net effect is small 
uncertainties (<TA g e ~ 0.1 dex) for ages less than 100 Myr, 
considerable uncertainties in the range 100 Myr to 1 Gyr 
(« 0.4 dex), and again quite small uncertainties for ages 
greater than 1 Gyr (w 0.2 dex). These numbers are similar 
for any evolutionary sequence with Z staxa < 0.004. For 
higher metallicities the age of many objects decreases, and 
the region of higher uncertainties at intermediate ages is 
shifted to smaller ages. Our conlusion is that the overall 
age pattern is quite insensitive to observational errors. 



5. Results 

We find an overall temporal age pattern that is quite ro- 
bust to changes in all model parameters except Z s tars and 
the lower mass limit. The latter matter only for ages com- 
parable to the lifetime of stars with mass equal to Mi QW , 
and M\ ow = 2M@ is unrealistic for ages > 1 Gyr. Despite 
some claims of IMFs devoid of low mass stars in starburst 
regions, there is no direct evidence for such, and we mainly 
consider such evolutionary sequences for comparison. 



In Fig. [7| we show age and mass histograms for a se- 
lection of evolutionary sequences with E(B — V) as a free 
parameter and fixed E(B — V) = 0.05. In accordance with 
OBR98, there is evidence both for young (less than 100 
Myr) intermediate (a few Gyr) and old (more than 10 
Gyr) clusters if a low metallicity is adopted. If a high 
metallicity is adopted, the net effect is to compress the 
age span, but retaining the main features. A metal-rich 
stellar component provides a poor fit to the bluest ob- 
jects which must therefore have a metallicity close to that 
derived for the nebular gas. Adopting a nebular metallic- 
ity of Z gas = 0.002, provides for a significantly better fit 
for young clusters, which give us confidence in our way of 
handling the nebular emission in Z2001. 

The shape of the mass histograms are model dependent 
(due to the different mass normalisations inherent to each 
IMF), but for most parameter choices we find that the 
maximum mass is higher for clusters older than 1 Gyr, and 
more mass seems to be found in relatively old than young 
clusters. All evolutionary sequences give a mass function 
(MF) with many clusters with M > 10 6 M Q , showing that 
this galaxy is rich in massive globular clusters. 

Comparing the free fit E[B - V) with the E(B - V) = 
0.05 results (top two panels in Fig. 01, one sees that peaks 
in the age distribution tend to be displaced to lower ages 
when E(B — V) is allowed to vary freely. This effect is 
notable for the peak at an age of a few Gyr. The peak is 
there for all metal-poor evolutionary sequences but occur 
at 0.2 dex lower age for the free fit E(B — V) cases. At 
these ages, few sources have accurate photometry in more 
than 2 bands (v and i) and the fixed E(B — V) results 
are probably more robust and less prone to systematic 
errors. However, the overall results in terms of age and 
mass distributions are insensitive to the adopted extinc- 
tion treatment, but this is so only because we have a priori 
information that the extinction must be low. 

When adopting a fixed E(B — V) = 0.05 (as in 
OBR98) we find a good overall agreement with the result 
in OBR98, except for ages below 10 Myr where the evolu- 
tionary sequences used in OBR98 had no time resolution. 
We do not find the peaks at 100 and 600 Myr reported by 
OBR98, although we do find objects with these ages. 

The age pattern is sensitive to the adopted star forma- 
tion time-scale tsf at ages comparable to t$f (compare 
top-left and mid- left panels in Fig. . The IMF influences 
the ages, but the effect is small at high ages; the effect on 
masses is more severe (upper and middle right panels in 
Fig.0. Both the a = 1.35 and 2.85 evolutionary sequences 
have problems at high ages since the predicted mass to 
light ratios (for 15 Gyr: M/L v = 47 and 8.7, respec- 
tively, as measured by the total consumed mass, whereas 
a = 2.35 give M/L v = 6 2 or M/L v = 4.2 if count- 
ing only stars and remnants) is much higher than those 
observed in old Galactic globular clusters (M/Ly ~ 2, 
Ashman & Zepf 1998). Gas parameters are only significant 
for ages smaller than or comparable to a few times tsf- 

2 A Scalo98 IMF would give a factor of ~2 lowevM/L v 
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Fig. 7. Comparison of output age and mass distributions for different model assumptions. In the left panels, the 
extinction was allowed to vary freely in the range E(B — V) S [0,0.25], and the resulting exinction distribution is 
indicated by the inset in the mass histograms. In the right panels a fixed extinction of E(B — V) = 0.05 was used. 
All evolutionary sequences, except for the lower right one, was calculated with Z2001 code, and have the following 
parameteers in common: Mass range M — 0.08 to 120M Q , Z gas = 0.002, and covering factor 1.0. Other parameters: 
top left: Salpeter IMF (a — 2.35), Z stars = 0.001 and instantaneous burst, middle left: a = 2.35, Z stais = 0.001 and 
constant SFR for 10 Myr. lower left row: a = 2.35, Z sta rs = 0.020 and constant SFR for 10 Myr. top right: Salpeter 
IMF (a = 2.35), Z stars = 0.001 and instantaneous burst, middle right: a — 1.35, Z staTS = 0.001 and instantaneous 
burst, lower right: PEGASE.2 model a = 2.35, Z sta , rs = 0.002 and instantaneous burst. 



The found E(B — V) distribution is relatively insensitive 
to the adopted parameters. Within each parameter con- 
figuration we have analysed the correlation between the 
fitted age, mass, residual and E{B — V). Mass and age are 
corellated (as expected since M/L is a strong function of 
age), but not the other quantities. 

We have performed correlation analysis on the out- 
put from the different parameter configurations. All evo- 
lutionary scqucncccs with a low stellar metallicity and 
A/iow = 0.08M Q show a high degree of mutual correla- 
tion, corroborating the results in Fig. The correlation 



with more metal rich sequences is significantly weaker but 
these parameter configurations viewed as a group show a 
high degree of internal correlation. Masses, RMS residu- 
als and E(B — V) values between different model runs are 
tightly correlated whenever the age-correlation is high. 

5.1. Comparison between Z2001 and PEGASE.2 

Figure also show the results from one of the PEGASE.2 
parameter configurations. We have already noted that the 
different treatments of the nebular component result in 
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very different behaviour during the first 10 Myr, and that 
PEGASE.2 cannot account for the most extreme colours. 
Also at later ages there are some differences between the 
two models. For a fixed E(B — V) the qualitative agree- 
ment for ages larger than 50 Myr is fair (Figure 0) , al- 
though PEGASE.2 give lower ages on average. For the 
E(B — V) free parameter cases, the pronounced peaks at 
~ 1 and 15 Gyr are replaced by a broader distribution 
for the PEGASE.2 model. Comparing the tracks in Fig. 
and 0] it is not surprising that we get different results, 
as the PEGASE.2 model show loops in the bvi diagram 
(Fig.[SJ) between 10 and 500 Myr and a redder colour at 1 
Gyr. For both Z2001 and PEGASE.2, a potential danger 
with a free fit E(B — V) is that features in the evolution- 
ary sequences, e.g. loops and kinks, may act as attractors. 
The degree of correlation between results from Z2001 and 
PEGASE.2 with similar parameter settings is high for age 
and mass but smaller for E(B — V) and RMS residuals. 

6. Resulting constraints on the starburst physics 

In agreement with the spectroscopic analysis, we find that 
young objects (age < 20 Myr) are better reproduced 
by metal-poor (Z stais < 0.004) stellar components. At 
somewhat higher ages, a slightly higher stellar metallic- 
ity (Z stais = 0.008) is also allowed and in some cases gives 
a slightly better fit. For ages on the order of 100 Myr 
and above, results are inconclusive regarding metallicity. 
However notable, if a higher metallicity is adopted, much 
fewer objects with ages in excess of 1 Gyr are found. 

We next turn our attention to the upper mass limit 
and slope of the IMF, which is important for regulating 
the ionising flux during the early stages of a burst. Hence 
it influences the UV-continuum and the ionised gas spec- 
trum. The combination of red b— v, and blue v—i colours in 
the colour-colour diagram at low ages in many evolution- 
ary scenarios is a direct consequence of significant contri- 
bution from gas emission to the ^-filter (F555W) whereas 
the b and i filters (F439W and F814W) are free of strong- 
lines, although nebular continuum is present in all filters. 
All scenarios with an upper mass limit of M up = 20M Q 
fail to reproduce the objects with red b — v and blue v — i 
colours (Fig. |Sj . Can we try to pin down the minimum al- 
lowed upper mass limit? A necessary requirement must be 
that such a scenario can reproduce the cloud of objects in 
the upper left part of Fig. |SJ To investigate this issue, we 
ran a new set of evolutionary sequences with M up = 40, 
60, and 80M© using Z001. 

Let us first look at the Salpeter IMF. For a covering 
factor of 0.5 the model cannot reproduce the observations 
even if M up = 120M©. For a covering factor of 1.0, there 
is no problem and an upper mass limit of 80Af Q can be 
accepted (Fig.|SJ|. However, if we increase the filling factor, 
we can mimic the effect of increasing M up . Hence, even if 
a filling factor of unity appears unrealistic according to 
the current knowledge of the ISM (Kassim et al. Q.989), 
we cannot rule out that it is higher than our generally 
assumed value of 0.1. Our more conservative conclusion is 



therefore that for a Salpeter IMF, M up > 60M Q . This is 
in agreement with Bergvall (1985) who from the observed 
[Oiii]a5oo7/H/3 ratio found that A/ up > 5OM . 

For a = 1.35, a covering factor of 1.0 give very extreme 
colours (Fig. EJ, suggesting that cluster formation has al- 
ready ceased or that M up < 60M©. Even a covering factor 
of 0.5 give quite extreme colours at the earliest times, but 
in general offset (to bluer v — i colours) from the cloud of 
points. This offset could possibly be explained by an av- 
erage extinction of E(B — V) = 0.15, which is much more 
than the estimates presented in Sect 4.1. If the sources 
are allowed to have a range of covering factors, an upper 
mass limit as low as M up > 40M© would provide a decent 
fit and the extinction may be lower. If we adopt a — 1.35 
we must assume that either: i) there are plenty of objects 
with age 10 Myr, but very few ones younger than this, or 
ii) AI up < 60M Q , or iii) the covering factor is 0.5 and the 
extinction is higher than previous estimates. 

For a — 2.85, a covering factor of 1.0 in combination 
with M up = 120Mq is needed to at least touch upon the 
the cloud of young clusters (Fig. . Hence we would be 
forced to adopt the lowest possible ages together with uni- 
form high filling and covering factors. We regard this as 
unrealistic, and hence a — 2.85 is not our preferred IMF. 

Could the IMF vary from cluster to cluster? The issue 
whether the IMF is universal or not has been subject to a 
many discussions (see Kroupa IZ002| . Variations could ei- 
ther be due to environmental or stochastical effects. The 
environmental effects should be similar for all clusters in 
ESO 338-04, but in low mass clusters the upper mass range 
of the IMF may be badly populated and subject to sta- 
tistical noise. For a Salpeter IMF with 0.08-120 M , the 
optical/UV emission from young clusters may be affected 
by stochastical effects for cluster masses below 10 5 M© 
(Langon and Mouhcinc 2000), which is close to our me- 
dian mass for your young clusters. Hence we cannot ex- 
clude stochastical effects among the fainter young clusters 
which could mimic a variable IMF slope. But in this case 
the average must be close to a — 2.35. 

7. Clustered vs diffuse emission 

At first sight (Fig. U [2J it appears that the star burst re- 
gion is dominated by emission from compact star clusters. 
However from Table 3 it is clear that even the central 1 
kpc is dominated by diffuse emission. Of course, the dif- 
fuse emission is also likely to be clustered at some level 
and include open clusters, OB associations etc. 

The integrated luminosity distribution of clusters in 
ESO 338-04 appears to have a Gaussian shape with a peak 
at My ~ — 10 (OBR98) but this is primarily a conse- 
quence of the photometric completeness limit and a mix 
of ages with different M/L. In each age bin containing 
a significant number of clusters, the luminosity function 
(LF) is best described by a power law dN/dL cx L~@, with 
power law index (3 = 1.8 ± 0.15 down to the completeness 
limit. This is valid also for the F218W data. From this we 
can estimate, roughly, how large a fraction of the clusters 
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Fig. 8. Colour-colour diagram (v — i vs b — v) showing, 
as dots, the location of clusters with uncertainties smaller 
than or equal to the rectangle in the lower left. The me- 
dian uncertainty is indicated by the cross. The arrow in- 
dicates how objects would be repositioned after correction 
for a hypothetical extinction of E(B — V) = 0.25. Lines 
represent the evolutionary tracks of the models from 0.5 
Myr (circles), 10 Myr (triangle), 1 Gyr (square) to 15 
Gyr (Star). The predictions at the youngest age are la- 
belled with the assumed upper mass limit. Parameters 
common for all scenarios: instantaneous burst, covering 
factor = 1.0, a = 2.35, Mi ow = 0.08M o , Z gas = 0.002, 
Z s tars = 0.001. 
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Fig. 9. Colour-colour diagram (v — i vs b — v) showing the 
location of clusters and evolutionary scenarios with dif- 
ferent IMF slope. Symbols have the same meaning as in 
Fig- [HI The scenarios are labelled with their assumed IMF 
slope: a — 1.35, 2.35 and 2.85. Other parameters common 
for all scenarios: instantaneous burst Afi ow = 0.08M Q , 
M up = 12OM , Z gas = 0.002, Z stars = 0.001, covering 
factor = 1.0. 



that fall below the completeness limit. We estimate, that 
if we correct for clusters below the detection threshold, 
the clusters would make up 40% of the light in F218W. 
However, although the number of young clusters is small 
outside the central kpc, diffuse Ha emission is present out 
to much greater radius (Bergvall & Ostlin 2002). Is this 
gas ionised in situ by a more diffuse population of young 
stars, or by leakage of Lyman continuum photons from 
the central burst? The Ha emission inside r = 1.075 kpc 
corresponds to 80% of the total Ha emission in ESO 338- 
04, which may bring the fraction of the total UV emission 
from the galaxy associated with young clusters, back down 
close to 30 % if the ionisation is in situ. Hence, the frac- 
tion of 2000A emission, and hence young stars, coming 
from young clusters should be in the range 30 to 40%. 

Meurer et al. l|TM5|) found only 16% of the 2200 A 
emission in ESO 338-IG04 to come from clusters, but this 
was estimated in a different way (using de-convolution 
of the pre-refurbishment HST/FOC images). The aver- 
age value for nine starbursts in Meurer et al. (199SJ) is 
20%. We have simply integrated the emission coming from 
clusters (including a significant aperture correction) in all 
passbands and compared to the integrated emission over 
a circularly defined starburst area. Also our integrated 
starburst magnitude (see Table 0J) is brighter by ^1 mag- 
nitude compared to Meurer et al. Ijl995|) . Apparently, our 
WFPC2 observations have been more sensitive and the 
fact that we get similar numbers for the other passbands, 
and find a monotonic increase in the cluster fraction with 
shorter wavelength, give us full confidence in our numbers. 

Zepf et al. |2001)) found that star clusters are respon- 
sible for 19% of the blue light in the the gas-rich merger 
NGC 3256, similar to our value of 21%, but most other 
mergers have a smaller fraction of their optical light com- 
ing from young clusters (Ostlin 2000), e.g. in the Antennae 
it is less than 10% (Whitmore & Zhang EES)- 

7.1. The ionisation energy budget 

,/From the Z2001 evolutionary scenarios, we have calcu- 
lated for each cluster, the predicted output of Lyman con- 
tinuum photons, and hence the predicted total Ha output 
associated with young clusters. The values we get have 
been scaled up with a factor ~ 1.3 to account for clusters 
fainter than the detection threshold. This correction has 
been applied to the estimates presented below. 

The total observed Ha flux from ESO 338-04 is 
log(Ha) = 41.7 ± 0.2 (Ostlin et al. 1999) which has only 
been corrected for Galactic extinction, but since the inter- 
nal extinction is estimated to be very small, E(B — V) < 
0.05, it will serve our purpose as an estimate. If the clus- 
ters make up 30 to 40% of the present star formation, they 
should then produce log(Ha) = 41.2 to 41.3(±0.2). 

A standard Salpeter metal-poor scenario with instant 
burst produces log(Ha) = 41.1. Hence some emission is 
lacking, but it is within the uncertainty. A constant SFR 
for 10 Myr gives 42.1, i.e. much larger than the observed 
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value. With a constant SFR for 5 Myr we get 41.9 and with 
an exponentially decaying SFR (e-folding time 3 Myr) we 
get 41.6. An IMF with a = 1.35 give similar values to a = 
2.35. Hence these results support that the star formation 
time-scale is short, a few Myr at most. 

8. Star formation history as derived from clusters 

The formation of a massive star cluster is by necessity as- 
sociated with star formation in the galaxy. In the present 
starburst in ESO 338-04, 30 to 40% of the newly formed 
stars are found in clusters. Are compact/massive star clus- 
ters formed as a natural consequence of star formation, 
or are special conditions, e.g. the high ISM turbulence 
prevailing in merging galaxies or circum- nuclear rings, re- 
quired? In the former scenario, clusters form now and then 
as long as star formation is active. Star formation is nor- 
mally a collective phenomenon, but only rarely (e.g. the 
Arches cluster, Serabyn et al.[T22H| do we find young stars 
to form self-gravitating systems in our Galaxy. This could 
be a consequence of a low star formation rate density 

The ability of young clusters to develop into globu- 
lar cluster-like systems is connected to their environment: 
Whether the clusters will be gravitationally bound or not, 
depends on e.g. the ISM density, pressure and metallic- 
ity Moreover, the ability for clusters to survive disruption 
depends on their environment (Meylan & Heggie 1997). 

For a given cluster age, the photometric completeness 
limit corresponds to a given stellar mass, producing a time 
dependent mass detection limit. The fraction of clusters 
in each age bin brighter than the limiting magnitude de- 
pends the cluster LF, a quantity which has been subject 
to some debate during the last years. The shape of the 
LF of old GCs in our Galaxy is Gaussian with a peak of 
log(M/M ) = 5.3 and dispersion 0.5 (Fall & ZhangEIIffl) • 
On the other hand, the LF and MF of young clusters, e.g. 
in the Antennae (Whitmore et al. 1999) and ESO 338-04, 
is described by a power law without any sign of a turnover 
at a characteristic mass ~ 10 5 M Q . As clusters age and 
are subject to internal and external dynamical effects, the 
lower part of the mass spectrum will be de-populated and 
the MF converted from power law to Gaussian (Fall and 
Zhang, 2001). However, the initial MF may still be a power 
law. Whether clusters survive or not, the power-law MF 
measures the mass initially formed in star clusters. 

If the initial cluster mass function is universal and 
known, it is straightforward to estimate what fraction of 
the clusters at a given age that fall below the detection 
threshold. However, in practise the mass integral for a 
power law with unbounded mass intervals diverges, and it 
is neccessary to fix the mass range, just like for the stellar 
IMF. Unfortunately, these quantities are not observation- 
ally known. One way would be to self consistently solve 
for the upper mass limit, but this only works when the 
MF is well populated. Moreover, this violates the assump- 
tion of a universal cluster MF and the same applies if the 
upper mass limit is taken as the observed maximum mass 
in the age bin in question. Our finally adopted method 



is quite crude and assumes a universal power-law initial 
mass function with a mass range 10 3 Mq to 5 x 1O 6 M . 

The time dependent correction factor for clusters be- 
low the detection treshold can now be calculated and 
ranges from 1.25 at low ages, where we see most of the 
mass (at least in systems with small reddening), to ~ 5 
at 15 Gyr where we see only about 20%. Assuming a 
Gaussian initial cluster mass function (e.g. to simulate 
that many low mass clusters get disrupted) would give 
similar numbers. It is not vitally important if clusters are 
gravitationally bound or not. Our estimate will always 
give a lower limit to the cluster SFR. 

For each scenario, we have included the uncertainty of 
±0.15 in the power law index f3 of the LF and an uncer- 
tainty in the completeness limit of 0.25 mag. The com- 
pleteness limit used is ~ lmag brighter than the photo- 
metric limit. This is to get closer to the real 100% com- 
pleteness, and also removes the most uncertain objects. 

8.1. The young starburst 

Figure [TOI shows the SFR vs time during the last 100 Myr 
as derived from clusters for 4 evolutionary scenarios. The 
patterns look close to identical. For instantaneous burts 
the main features are: a small peak at < t < 3 Myr, 
followed by a decrease for 3 < t < 6 Myr, after which the 
main peak occur in the interval 6 < t < 10 Myr. Another 
significant peak occurs at 20 < t < 30 after which the 
cluster SFR drops slowly but remains positive until 40 
Myr. 

For longer bursts (constant SFR during 10 Myr, 3rd 
panel in Fig. llOJI the highest peak is shifted to the interval 
10 < t < 15 Myr, while the peak at 25 Myr is not affected 
(see third panel in Fig. I10[) . Changing the IMF does not 
affect the pattern, only the scaling of the SFR. 

The relatively low SFR in the first bin (0 < t < 3 Myr) 
is not a consequence of a few number of clusters, but is due 
to their low estimated total mass. The scarcity of objects 
with ages near 5 Myr can be seen also in Figs. |H1 and [5j 
many objects are clustered around v — i = 0, 6 — v = 
and v — i = — 1,6 — v = 1, but few in between. Clearly, 
at such small ages with rapid stellar evolution, ages will 
be sensitive to details in the stellar models. The dip at 15 
Myr could possibly be attributed to features in the multi- 
dimensional colour plane described by the scenarios. 

Whether or not the SFR went down to zero some 15 
Myr ago, the star formation activity began to rise already 
40-50 Myr ago, and appears to have peaked a little less 
than 10 Myr ago. Hence the starburst in ESO338-IG04 is 
40-50 Myr old and still active. This estimate is based on 
clusters selected from their bright optical and UV emis- 
sion. It is possible that a fraction of the youngest sources 
are hidden in dusty cocoons, that should be visible as mid- 
IR sources like those found in the BCG He2-10 (Vacca et 
al. 12002(1 . but see discussion below. 

If there is a time-lag of a few Myr in the formation 
of the most massive stars (see discussion above), some of 
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Fig. 10. The star formation rate derived from clusters as 
a function of time during the past 100 Myr for 4 different 
evolutionary scenarios, all with Z gas = 0.002, Z s t a rs = 
0.001 and mass range 0.08 - 120 M@. First (top) panel: 
Salpeter IMF (a = 2.35) with instantaneous burst and 
covering factor = 1.0. Second panel: Salpeter IMF with 
instantaneous burst when the covering factor is allowed to 
vary. Third panel: Salpeter IMF with a burst duration 
of 10 Myr and covering factor 1.0. Fourth panel: a = 
1.35 IMF with instantaneous burst and covering factor 
0.5. Error bars in the X-direction shows the width of the 
time window, and points are connected with dotted lines 
to guide the eye. Error bars in the Y-direction account for 
the uncertainty of ±0.15 in the power law index (3 of the 
LF and an uncertainty in the completeness limit of 0.25 
mag. 

the clusters in the 5-10 Myr peak would instead be very 
young (~ 1 Myr) and the difference in SFR between the 
0-5 and 5-10 Myr intervals diminished. 

8.2. The SFR viewed over a cosmological time scale 

Judging from the scenarios with a metal-poor stellar com- 
ponent, there is no evidence for a significantly enhanced 
SFR 50-1000 Myr ago, as derived from the observed clus- 
ters. However, if we set -2T s tars > 0.008, the cluster SFR 
remains high over the past Gyr although all the peaks at 
t < 100 Myr remain. What happens is that many clusters 
are found to be younger if the metallicity is increased. 
Figure shows the cluster SFR over the past 15 Gyr 
for three different metallicities (Z star s = 0.001, 0.008 and 



0.020). In addition, we show the result from a PEGASE.2 
evolutionary scenario. The present starburst is still ev- 
ident in the first bin (0 < t < 500 Myr) but has de- 
creased in amplitude since the SFR has been averaged 
over a longer time-scale. For the Z staTS = 0.001 scenario, a 
prominent peak is found at 1.5 Gyr, and significant cluster 
formation occurred also between 2 and 5 Gyr ago, which 
was found also by OBR98. Assuming that all GCs older 
than 10 Gyr formed during a short time scale (r < 1 Gyr) 
would produce a peak comparable to that found at 1.5 
Gyr. The event 1.5 Myrs ago must have been an impres- 
sive starburst. 



Bergvall & Ostlin (2002) reported that the underly- 
ing host galaxy may be significantly more metal-rich than 
the nebular gas. In this case one might suspect that older 
clusters could be metal-rich. For higher metallicities, the 
1.5 Gyr feature remain but decrease in amplitude. The 
14 Gyr peak, on the other hand, becomes stronger. This 
is due to the higher M/L at high ages in these scenar- 
ios, making both the mass detection limit and the masses 
of individual clusters higher. The unknown metallicity of 
older clusters and the small number of objects in each age 
bin presents the most important limitation for estimating 
the past SFR from clusters. 



8.3. The cluster mass vs total galaxy mass 

Summing up the total cluster SFR over 15 Gyrs we get 
(SFR) = 0.02Af Q /yr, or a total mass formed in clusters 
of ~ 3 x 1O 8 M0. The total stellar mass of this galaxy is 
estimated to be ~ 4 x 10 9 M Q (Ostlin et al. lMiTl Bergvall 
& Ostlin 2002). Hence the clusters would be responsible 
for almost 10% of the total star formation. This is orders 
of magnitude higher than for most galaxies, for which the 
typical value is 0.3% (McLaughlin IT333I I2TTO . For in- 
stance, the total mass of the Milky Way globular cluster 
system is estimated to be on the order of ~ 6 x 1O 7 M0. 
As noted in Sect. 4.1, if the actual IMF is more similar 
to Scalo98 than Salpeter, we will overestimate the clus- 
ter masses. In addition, we count the total consumed gas 
mass (including the gas returned from stellar winds and 
explosions), but these two effects combined should not be 
greater than a factor of 2. Old clusters may have suffered 
loss of low mass stars through dynamical evolution, but 
also this effect should be not greater than a factor of 2 
in view of the mass to light ratios for Galactic globular 
clusters (Ashman &; Zenf I199%| . Taking also the mass of 
gas M(Hl) = 3 x 10 9 Mq into account we are still an 
order of magnitude above 0.3%. If we count only the ac- 
tually detected clusters we still get a total cluster mass of 
9 x 10 7 A/ Q for Z stars = 0.001 (> 7 x 10 7 for any Z stars ), 
or about one third of the total implied cluster mass. This 
would still give a mass fraction ^> 0.3% if accounting for 
the possible systematic effects mentioned above. 

The conclusion that this galaxy is, and has been over 
a cosmological time-scale, a very efficient cluster former is 
unavoidable. The open question is why? Tashiro and Nishi 
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Fig. 11. The star formation rate in clusters as a function 
of time during the past 15 Gyr for three Z2001 evolution- 
ary sequences with different Z s tars (0.001, 0.008, and 0.02), 
and one PEGASE. 2 sequence with Z atars = 0.002. All se- 
quences have a = 2.35, mass range 0.08 - 120 M Q and 
Z gas = 0.002. The top panel shows the correction factor 
due to clusters below the detection limit for the different 
metallicities. 



( 2000) suggest that the formation of gravitationally bound 
clusters may be favoured in low metallicity environments. 

Bergvall & Ostlin (2002} found that the stellar metal- 
licity of the host galaxy is higher than that of the ISM, 
which could be explained if a metal-poor gas-rich dwarf 
galaxy or an H I cloud has been accreted. They speculate 
that the unusually high metallicity, in view of the low lu- 
minosity, arise because the galaxy is rich in dark matter 
(see Sect 19.1( 1. which has enabled it to retain all its metals. 
If this is true, the total mass fraction (including dark mat- 
ter) of globular and akin clusters may not be higher than 
in other galaxies. Another possibility if one postulates a 
universal cluster formation efficiency is that for some rea- 
son the survival chances are greater in ESO 338-IG04, but 
then this should be a general trend among dwarf galaxies. 



9. Spatial propagation of the starburst 

The distribution of young clusters can be used to investi- 
gate if propagating star formation have occurred. In Fig. 
I12l we show the spatial age pattern for young clusters for a 
standard metal-poor Salpeter scenario (our favoured one 
for young objects) with a varying covering factor. There 
are indeed some striking patterns present in this image. 

The youngest objects (< 3 Myr), coded in cyan, tend 
to prefer a mix of covering factors. They are preferentially 
found in a ring like structure centred on X=270, Y=305 
(capital X and Y refer to Planetary Camera pixel values, 
each pixel has a size of 8.3 pc, see Fig. ll2f) and in a group 
near X=320, Y=400. The former structure have in general 
a best fitting covering factor of 1.0, and the latter 0.5. 

Slightly older objects (3 < Age < 8 Myr, coded in 
blue), are found in the latter group and also near the star- 
burst centre (X=290, Y=340). Objects with ages 8 to 30 
Myr are found in the centre and in a group near (X=300, 
Y=260). Older objects (30-100) Myr are not strongly clus- 
tered, but are preferentially found in the upper half. 

The ring- like structure of very young (< 3 Myr) ob- 
jects have a major axis diameter of 0.8 kpc. A velocity of 1 
km/s correspond to 1 pc/Myr. The Ha velocity dispersion 
is <7Ha = 43 km/s (Ostlin et al. 2001) implying a turbulent 
half ring crossing time of 10 Myr. Hence the coeval for- 
mation of this ring can not have been triggered by a wave 
travelling at the ISM turbulent velocity, but perhaps by a 
supernovae driven wind. The grouping of very young ob- 
jects in the upper part is more compact and consistent 
with a coeval formation if random motions of 0.1 kpc are 
considered. To the right of this grouping several diffuse 
emission-line regions without compact clusters are present 
(see Fig. ^ , which may possibly be the young sites of star 
formation. The ~ 10 Myr (yellow) grouping in the lower 
part (X=300, Y=260) has a similar dimension. Their age 
is consistent with a localised formation followed by diffu- 
sion by random motions. In this region we find an inter- 
esting boundary between ~ 10 and ~ 25 Myr old objects, 
although one can not exclude that a dust filament is re- 
sponsible for the apparent ridge of ~ 25 Myr objects. The 
centre (X=290, Y=340) is close to the centre of the outer 
isophotes and could therefore represent the central poten- 
tial of the galaxy. Here cluster formation appears to have 
been continuing over the last 40 Myr. The present star- 
burst may have started in the region slightly west (up in 
Fig.ll2|) of the centre where almost all of the sources older 
than 30 Myr are located. A possible scenario is that the 
upper and lower 8 to 30 Myr groups represent approaching 
complexes that have triggered the young shell-like struc- 
ture. Indeed, the Ha velocity field show that the eastern 
part of the galaxy (down-left in Fig. ll2[) may be falling in 
towards the centre (Ostlin et al. 2001). Moreover, the Ha 
velocity dispersion is higher east of this region (cthq = 70 
km/s) and there are suggestions of a perpendicular ve- 
locity component. Further out towards east (outside the 
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Fig. 12. Spatial distribution of young clusters using best fitting E(B — V) and Salpeter IMF with mass range 0.08 
to 120 Mq, and metallicity Z s tars = 0.001 and Z gas = 0.002. The covering factor was allowed to vary between 0.5 (+ 
signs) and 1.0 (dots). The size of each symbol is proportional to the logarithm of the mass of a cluster and the colour 
depends on the age. The x- and y- coordinates correspond to pixel numbers in the Planetary Camera and the region 
showed is 14.6" x 14.6", or 2.6 x 2.6 kpc. 



Planetary Camera FOV shown in Fig. ^| we find some 
diffuse emission-line regions and clusters 3 . 

Although projection effects and the velocities of indi- 
vidual clusters are not known, our results provide evidence 
for the spatial propagation of the starburst over the last 
30 Myr. Supersonic triggering appears to be necessary to 
explain some of the youngest features. 

The star formation pattern and the conclusions remain 
unaffected if a fixed E{B — V) — 0.0 or 0.05, or a fixed 
covering factor, is adopted. The same applies to changing 
the IMF: ages would change slightly, but not the pattern. 
If massive star formation would be associated with a time- 
lag, as discussed above, most features would remain, e.g. 



3 Some blueish clusters in this tail ( Ostlin et al. 2001, 
Bergvall & Ostlin 2002) are visible in the Wide Field chip 3, 
but their photometry is too uncertain to assess their ages. 



the young ring, but it would be a few Myrs older, but 
still too young for a sound wave to have triggered their 
formation. Some of the ~10 Myr objects (in yellow) would 
now rather be just a few Myr. 

9.1. Spatial distribution of old globular clusters 

For old (> 0.8 Gyr) clusters the time-scale is too long for 
the position of individual clusters to have any correspon- 
dence with their place of birth. Instead, we examine the 
spatial distribution of the ensemble of old clusters (Fig. 
1 1 31) . Intermediate age (0.8 < age < 5 Gyr) clusters (open 
shaded symbols in Fig. I13|) have a flattened distribution 
with minor to major axis ratio b/a ~ 0.5, with the same 
position angle as the galaxy as a whole. The isophotes 
of ESO 338-04 have an axis ratio b/a = 0.48 (Bergvall 
& Ostlin |2TJ0"2|) . The colours of the host galaxy at radii 
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Fig. 13. Spatial distribution of old (age > 5 Gyr, filled 
black symbols) and intermediate age (0.8 < age < 5 Gyr, 
open shaded symbols) clusters with size proportional to 
mass. The small dots younger (< 0.8 Gyr) objects. The 
axes show the spatial coordinates, in units of PC pixels. 
The region delimited by a dashed line (X < 50, Y < 50) 
lies at the border to the wide field camera CCDs and is 
subject to vignetting, and therefore not included in our 
analysis. 



sufficiently large that the starburst contribution is negli- 
gible, indicate an age in the range 2-8 Gyr and a relatively 
high stellar metallicity. That the intermediate age globu- 
lar cluster system has an identical flattening suggests a 
common origin or, if the lower age of the clusters is taken 
at face value, that these formed in the existing gravita- 
tional potential of the host galaxy. If a higher metallicity 
is adopted for the clusters, the age difference would in- 
crease. 

On the contrary, clusters older than 5 Gyr (filled 
black symbols in Fig. 113(1 show no apparent flattening 
(b/a w 1.0) suggesting the presence of an old spherical 
halo, and that they are indeed metal-poor. The number of 
objects is too small to study the radial distribution in de- 
tail, but the effective radius (containing half the number of 
>5 Gyr clusters) is ~ 1.4 kpc. If a Gaussian LF similar to 
the one for Galactic GCs (MP eak = -7.2, a = 1.3 mag) is 
assumed, the total number (including non detected ones) 
of old halo clusters should be ~50. For a typical GC rich- 
ness of spherical systems (T = 5, see Ashman & Zepf[l998 
for a definition of T) 50 GCs would then imply a mass of 
the halo of 10 10 Af Q . This thought experiment supports 
the idea of a dark matter-rich halo discussed in Sect. 18.31 



10. Discussion - Star formation and dusty 
molecular clouds 

Star formation in our Galaxy and elsewhere is associated 
with dusty giant molecular clouds. Massive starbursts like 
M82 and Arp220 are reddened by several tens of magni- 
tudes of optical extinction (Genzel et al. H998|) . and a large 
fraction of the UV power is absorbed and re-radiated in 
the infrared. However, the extinction inferred from young 
star clusters in M82 is only of the order of Ay = 1 (de Grijs 
et al. I2001fl compared to the 50 magnitudes derived from 
IR data (Forster Schrcibcr et al. 2001). In the Antennae, 
the youngest clusters have high extinction Ay — 3 to 5 
(Zhang et al|2001), and are spatially correlated with the IR 
emission, whereas the nuclei have Ay ~ 50. Obviously, op- 
tical selection systematically picks out sources with lower 
than average extinction. 

The situation in ESO 338-04 is quite different. Looking 
at Fig. 21 and we see no tendency for the youngest ob- 
jects to be reddened. A few sources lie to the upper right 
compared to the locus of the evolutionary sequences at the 
lowest age, as if they were reddened by E(B — V) < 0.25. 
However, there are more sources to the down-left, and the 
observations are consistent with a dust free population 
plus random effects (integer number of massive stars, and 
small variations is gas filling and covering factors) plus 
photometric errors. Although our sample is biased towards 
optically bright clusters it is nevertheless surprising that 
we do not see any young sources just about to come out of 
their dusty cocoons in which they, according to common 
wisdom, should have formed. Note that there are nebu- 
lar emission dominated regions without associated point 
sources west of the centre, which may represent young ob- 
jects, still partly embedded. 

The lifetime of Galactic giant molecular clouds against 
photo- dissociation and ionisation is on the order of 30 
Myr once massive star formation commence (Williams and 
McKee 1997). If we insist that star formation is associ- 
ated with dusty molecular clouds, then the lack of young 
reddened sources suggest that in metal-poor starbursts 
with high ISM turbulence (Ostlin et al. 1200 l|l and high 
star formation efficiency, this time scale is much shorter. 
ESO 338-04 is a hot IRAS source {S 6 a/S W0 = 1.7) but 
the FIR luminosity is modest L\-r = 5.8 x 10 9 L Q and the 
IR to blue luminosity, Lir/Lb = 1.0 (solar units), is very 
low with respect to its /60//100 ratio (Bergvall et al.|2000 
Dulzin-Hacyan et al. 1988). The derived dust mass is only 
of the order of 1O 5 M (see also Calzetti et al lTMB) . This 
implies small amounts of dust heated to high tempera- 
tures, and fits well with the short time-scale scenario. 

A short time-scale for molecular cloud destruction 
would be consistent with the difficulty to find CO emis- 
sion in metal-poor dwarf starbursts (Kunth & Ostlin 
2000) since not only the CO to H2 conversion factor 
may be low, but the gas would stay in molecular form 
a shorter time. For ESO 338-04, L CO {i-o) < 7 x 10 27 W, 
so \og{Lcol 'Lib) < —8.4 (Bergvall, unpublished), imply- 
ing Mh 2 < 7 x 10 6 M Q if a standard CO to H 2 conversion 
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factor is used. As the mass of young (< 10 Myr) stars 
(in clusters and diffuse) is ~ 1O 7 M0, and star formation 
efficiency must be well below 100%, the molecular gas dis- 
sociation time scale should be a few Myr at most. 

While we see several very young (~1 Myr) clusters, 
Fig. 1 101 shows that the global star formation intensity may 
have passed its peak. This combined with the tendency 
to pick out sources with low extinction suggest that some 
star formation could be hidden after all. Of course we may 
happen to observe this galaxy just as the last reservoirs of 
dusty molecular gas are depleted, but we would still need 
to postulate a short dust destruction time scale. 

If the time scale for molecular cloud destruction is un- 
usually short, it is possible that a young cluster would still 
not have reached its peak UV luminosity when becoming 
naked. This could give us an additional selection effect to 
account for the lack of young reddened clusters. 

The cluster richness implies a high cluster and star 
formation efficiency (see Sect. 8.3), which could be a nat- 
ural cause for the short dust destruction time-scales since 
the UV radiation density in the ISM is high. It is surpris- 
ing that the star formation efficiency can be so high in a 
galaxy with a metal-poor ISM and little dust. 

11. Conclusions 

Using multicolour photometry from 2000 to 8000A, we 
have analysed the temporal and spatial evolution of star 
clusters and the starburst in the Luminous Blue Compact 
Galaxy ESO 338-04. Our main results are: 

- Using an extensive set of spectral evolutionary syn- 
thesis scenarios and Monte-Sarlo simulations we conclude 
that the overall age distribution found is very robust. 
Massive star clusters with ages from 1 Myr up to the max- 
imum age considered in the models (15 Gyr) are present 
in this galaxy. 

- A careful treatment of the nebular emission is crucial 
for modelling the properties of young star clusters. 

- In accordance with previous studies, we find the star- 
burst in ESO 338-04 to be almost unaffected by dust ex- 
tinction. There is no evidence for a large population of 
young sources about to come out of their dusty cocoons. 
We speculate that this can be explained by a short time- 
scale, < 1 Myr, for destruction of dusty molecular clouds. 

- The cluster formation history have been used to con- 
strain the overall star formation history in this galaxy. The 
present starburst began 40-50 Myr ago and is still active. 
There are indications that the SFR has varied during time, 
and that it has already passed its peak. 

- Young clusters are responsible for 30 to 40% of the 
star formation and the starburst luminosity output in the 
UV, which is a high number, even for starbursts. 

- Seen over a cosmological time, this galaxy has been 
a remarkably efficient cluster former. The reason for this 
remains an enigma. The intermediate age population has a 
flattened distribution with the same flattening as the host 
galaxy. The old clusters have a spherical distribution. 



- The observations are well reproduced by a Salpeter 
IMF (a = 2.35) with an upper mass limit M up > 60M Q . 
A steeper IMF (a = 2.85) cannot reproduce the data. 
A flatter IMF (a — 1.35) cannot be excluded but requires 
additional constraints to provide a good fit. An ISM metal- 
licity close to the spectroscopic value is strongly favoured. 
Our results suggest that the formation time-scale for an 
individual cluster is at maximum a few Myrs. 

- We find evidence for spatially correlated formation 
of groups of clusters, some requiring supersonic triggering 
velocities, e.g. from supernovae winds. 
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